Xyloglucan deficiency leads to a reduction in turgor pressure and changes in cell wall properties, affecting early seedling establishment.

Xyloglucan is believed to play a significant role in cell wall mechanics of dicot plants. Surprisingly, Arabidopsis plants defective in xyloglucan biosynthesis exhibit nearly normal growth and development. We investigated a mutant line, cslc-Δ5, lacking activity in all five Arabidopsis cellulose synthase like-C (CSLC) genes responsible for xyloglucan backbone biosynthesis. We observed that this xyloglucan-deficient line exhibited reduced cellulose crystallinity and increased pectin levels, suggesting the existence of feedback mechanisms that regulate wall composition to compensate for the absence of xyloglucan. These alterations in cell wall composition in the xyloglucan-absent plants were further linked to a decrease in cell wall elastic modulus and rupture stress, as observed through atomic force microscopy (AFM) and extensometer-based techniques. This raised questions about how plants with such modified cell wall properties can maintain normal growth. Our investigation revealed two key factors contributing to this phenomenon. First, measurements of turgor pressure, a primary driver of plant growth, revealed that cslc-Δ5 plants have reduced turgor, preventing the compromised walls from bursting while still allowing growth to occur. Second, we discovered the conservation of elastic asymmetry (ratio of axial to transverse wall elasticity) in the mutant, suggesting an additional mechanism contributing to the maintenance of normal growth. This novel feedback mechanism between cell wall composition and mechanical properties, coupled with turgor pressure regulation, plays a central role in the control of plant growth and is critical for seedling establishment in a mechanically challenging environment by affecting shoot emergence and root penetration.

Stiffness transitions in new walls post-cell division differ between Marchantia polymorpha gemmae and Arabidopsis thaliana leaves.

Plant morphogenesis is governed by the mechanics of the cell wall-a stiff and thin polymeric box that encloses the cells. The cell wall is a highly dynamic composite material. New cell walls are added during cell division. As the cells continue to grow, the properties of cell walls are modulated to undergo significant changes in shape and size without breakage. Spatial and temporal variations in cell wall mechanical properties have been observed. However, how they relate to cell division remains an outstanding question. Here, we combine time-lapse imaging with local mechanical measurements via atomic force microscopy to systematically map the cell wall's age and growth, with their stiffness. We make use of two systems, Marchantia polymorpha gemmae, and Arabidopsis thaliana leaves. We first characterize the growth and cell division of M. polymorpha gemmae. We then demonstrate that cell division in M. polymorpha gemmae results in the generation of a temporary stiffer and slower-growing new wall. In contrast, this transient phenomenon is absent in A. thaliana leaves. We provide evidence that this different temporal behavior has a direct impact on the local cell geometry via changes in the junction angle. These results are expected to pave the way for developing more realistic plant morphogenetic models and to advance the study into the impact of cell division on tissue growth.

Mean-field theory approach to three-dimensional nematic phase transitions in microtubules.

Microtubules are dynamic intracellular fibers that have been observed experimentally to undergo spontaneous self-alignment. We formulate a three-dimensional (3D) mean-field theory model to analyze the nematic phase transition of microtubules growing and interacting within a 3D space, then make a comparison with computational simulations. We identify a control parameter G_{eff} and predict a unique critical value G_{eff}=1.56 for which a phase transition can occur. Furthermore, we show both analytically and using simulations that this predicted critical value does not depend on the presence of zippering. The mean-field theory developed here provides an analytical estimate of microtubule patterning characteristics without running time-consuming simulations and is a step towards bridging scales from microtubule behavior to multicellular simulations.

Context-specific functions of transcription factors controlling plant development: From leaves to flowers.

Plant development is regulated by transcription factors that often act in more than one process and stage of development. Yet the molecular mechanisms that govern the functional diversity and specificity of these proteins remains far from understood. Flower development provides an ideal context to study these mechanisms since the development of distinct floral organs depends on similar but distinct combinations of transcriptional regulators. Recent work also highlights the importance of leaf polarity regulators as additional key factors in flower initiation, floral organ morphogenesis, and possibly floral organ positioning. A detailed understanding of how these factors work in combination will enable us to address outstanding questions in flower development including how distinct shapes and positions of floral organs are generated. Experimental approaches and computer-based modeling will be required to characterize gene-regulatory networks at the level of single cells.

High-Throughput 3D Phenotyping of Plant Shoot Apical Meristems From Tissue-Resolution Data.

Confocal imaging is a well-established method for investigating plant phenotypes on the tissue and organ level. However, many differences are difficult to assess by visual inspection and researchers rely extensively on ad hoc manual quantification techniques and qualitative assessment. Here we present a method for quantitatively phenotyping large samples of plant tissue morphologies using triangulated isosurfaces. We successfully demonstrate the applicability of the approach using confocal imaging of aerial organs in Arabidopsis thaliana. Automatic identification of flower primordia using the surface curvature as an indication of outgrowth allows for high-throughput quantification of divergence angles and further analysis of individual flowers. We demonstrate the throughput of our method by quantifying geometric features of 1065 flower primordia from 172 plants, comparing auxin transport mutants to wild type. Additionally, we find that a paraboloid provides a simple geometric parameterisation of the shoot inflorescence domain with few parameters. We utilise parameterisation methods to provide a computational comparison of the shoot apex defined by a fluorescent reporter of the central zone marker gene CLAVATA3 with the apex defined by the paraboloid. Finally, we analyse the impact of mutations which alter mechanical properties on inflorescence dome curvature and compare the results with auxin transport mutants. Our results suggest that region-specific expression domains of genes regulating cell wall biosynthesis and local auxin transport can be important in maintaining the wildtype tissue shape. Altogether, our results indicate a general approach to parameterise and quantify plant development in 3D, which is applicable also in cases where data resolution is limited, and cell segmentation not possible. This enables researchers to address fundamental questions of plant development by quantitative phenotyping with high throughput, consistency and reproducibility.

A 3D gene expression atlas of the floral meristem based on spatial reconstruction of single nucleus RNA sequencing data.

Cellular heterogeneity in growth and differentiation results in organ patterning. Single-cell transcriptomics allows characterization of gene expression heterogeneity in developing organs at unprecedented resolution. However, the original physical location of the cell is lost during this methodology. To recover the original location of cells in the developing organ is essential to link gene activity with cellular identity and function in plants. Here, we propose a method to reconstruct genome-wide gene expression patterns of individual cells in a 3D flower meristem by combining single-nuclei RNA-seq with microcopy-based 3D spatial reconstruction. By this, gene expression differences among meristematic domains giving rise to different tissue and organ types can be determined. As a proof of principle, the method is used to trace the initiation of vascular identity within the floral meristem. Our work demonstrates the power of spatially reconstructed single cell transcriptome atlases to understand plant morphogenesis. The floral meristem 3D gene expression atlas can be accessed at http://threed-flower-meristem.herokuapp.com .

KATANIN and CLASP function at different spatial scales to mediate microtubule response to mechanical stress in Arabidopsis cotyledons.

Mechanical stress influences cell- and tissue-scale processes across all kingdoms. It remains challenging to delineate how mechanical stress, originating at these different length scales, impacts cell and tissue form. We combine growth tracking of cells, quantitative image analysis, as well as molecular and mechanical perturbations to address this problem in pavement cells of Arabidopsis thaliana cotyledon tissue. We show that microtubule organization based on chemical signals and cell-shape-derived mechanical stress varies during early stages of pavement cell development and is mediated by the evolutionary conserved proteins, KATANIN and CLASP. However, we find that these proteins regulate microtubule organization in response to tissue-scale mechanical stress to different extents in the cotyledon epidermis. Our results further demonstrate that regulation of cotyledon form is uncoupled from the mechanical-stress-dependent control of pavement cell shape that relies on microtubule organization governed by subcellular mechanical stress.

Molecular mechanism of cytokinin-activated cell division in Arabidopsis.

Mitogens trigger cell division in animals. In plants, cytokinins, a group of phytohormones derived from adenine, stimulate cell proliferation. Cytokinin signaling is initiated by membrane-associated histidine kinase receptors and transduced through a phosphorelay system. We show that in the Arabidopsis shoot apical meristem (SAM), cytokinin regulates cell division by promoting nuclear shuttling of Myb-domain protein 3R4 (MYB3R4), a transcription factor that activates mitotic gene expression. Newly synthesized MYB3R4 protein resides predominantly in the cytoplasm. At the G2-to-M transition, rapid nuclear accumulation of MYB3R4-consistent with an associated transient peak in cytokinin concentration-feeds a positive feedback loop involving importins and initiates a transcriptional cascade that drives mitosis and cytokinesis. An engineered nuclear-restricted MYB3R4 mimics the cytokinin effects of enhanced cell proliferation and meristem growth.

A multiscale analysis of early flower development in Arabidopsis provides an integrated view of molecular regulation and growth control.

We have analyzed the link between the gene regulation and growth during the early stages of flower development in Arabidopsis. Starting from time-lapse images, we generated a 4D atlas of early flower development, including cell lineage, cellular growth rates, and the expression patterns of regulatory genes. This information was introduced in MorphoNet, a web-based platform. Using computational models, we found that the literature-based molecular network only explained a minority of the gene expression patterns. This was substantially improved by adding regulatory hypotheses for individual genes. Correlating growth with the combinatorial expression of multiple regulators led to a set of hypotheses for the action of individual genes in morphogenesis. This identified the central factor LEAFY as a potential regulator of heterogeneous growth, which was supported by quantifying growth patterns in a leafy mutant. By providing an integrated view, this atlas should represent a fundamental step toward mechanistic models of flower development.

Tissue folding at the organ-meristem boundary results in nuclear compression and chromatin compaction.

Artificial mechanical perturbations affect chromatin in animal cells in culture. Whether this is also relevant to growing tissues in living organisms remains debated. In plants, aerial organ emergence occurs through localized outgrowth at the periphery of the shoot apical meristem, which also contains a stem cell niche. Interestingly, organ outgrowth has been proposed to generate compression in the saddle-shaped organ-meristem boundary domain. Yet whether such growth-induced mechanical stress affects chromatin in plant tissues is unknown. Here, by imaging the nuclear envelope in vivo over time and quantifying nucleus deformation, we demonstrate the presence of active nuclear compression in that domain. We developed a quantitative pipeline amenable to identifying a subset of very deformed nuclei deep in the boundary and in which nuclei become gradually narrower and more elongated as the cell contracts transversely. In this domain, we find that the number of chromocenters is reduced, as shown by chromatin staining and labeling, and that the expression of linker histone H1.3 is induced. As further evidence of the role of forces on chromatin changes, artificial compression with a MicroVice could induce the ectopic expression of H1.3 in the rest of the meristem. Furthermore, while the methylation status of chromatin was correlated with nucleus deformation at the meristem boundary, such correlation was lost in the h1.3 mutant. Altogether, we reveal that organogenesis in plants generates compression that is able to have global effects on chromatin in individual cells.

Limits and Constraints on Mechanisms of Cell-Cycle Regulation Imposed by Cell Size-Homeostasis Measurements.

High-throughput imaging has led to an explosion of observations about cell-size homeostasis across the kingdoms of life. Among bacteria, "adder" behavior-in which a constant size increment appears to be added during each cell cycle-is ubiquitous, while various eukaryotes show other size-homeostasis behaviors. Since interactions between cell-cycle progression and growth ultimately determine such behaviors, we developed a general model of cell-cycle regulation. Our analyses reveal a range of scenarios that are plausible but fail to regulate cell size, indicating that mechanisms of cell-cycle regulation are stringently limited by size-control requirements, and possibly why certain cell-cycle features are strongly conserved. Cell-cycle features can play unintuitive roles in altering size-homeostasis behaviors: noisy regulator production can enhance adder behavior, while Whi5-like inhibitor dilutors respond sensitively to perturbations to G2/M control and noisy G1/S checkpoints. Our model thus provides holistic insights into the mechanistic implications of size-homeostasis experimental measurements.

Cytoskeletal organization in isolated plant cells under geometry control.

The cytoskeleton plays a key role in establishing robust cell shape. In animals, it is well established that cell shape can also influence cytoskeletal organization. Cytoskeletal proteins are well conserved between animal and plant kingdoms; nevertheless, because plant cells exhibit major structural differences to animal cells, the question arises whether the plant cytoskeleton also responds to geometrical cues. Recent numerical simulations predicted that a geometry-based rule is sufficient to explain the microtubule (MT) organization observed in cells. Due to their high flexural rigidity and persistence length of the order of a few millimeters, MTs are rigid over cellular dimensions and are thus expected to align along their long axis if constrained in specific geometries. This hypothesis remains to be tested in cellulo Here, we explore the relative contribution of geometry to the final organization of actin and MT cytoskeletons in single plant cells of Arabidopsis thaliana We show that the cytoskeleton aligns with the long axis of the cells. We find that actin organization relies on MTs but not the opposite. We develop a model of self-organizing MTs in three dimensions, which predicts the importance of MT severing, which we confirm experimentally. This work is a first step toward assessing quantitatively how cellular geometry contributes to the control of cytoskeletal organization in living plant cells.

It's about time: Analysing simplifying assumptions for modelling multi-step pathways in systems biology.

Thoughtful use of simplifying assumptions is crucial to make systems biology models tractable while still representative of the underlying biology. A useful simplification can elucidate the core dynamics of a system. A poorly chosen assumption can, however, either render a model too complicated for making conclusions or it can prevent an otherwise accurate model from describing experimentally observed dynamics. Here, we perform a computational investigation of sequential multi-step pathway models that contain fewer pathway steps than the system they are designed to emulate. We demonstrate when such models will fail to reproduce data and how detrimental truncation of a pathway leads to detectable signatures in model dynamics and its optimised parameters. An alternative assumption is suggested for simplifying such pathways. Rather than assuming a truncated number of pathway steps, we propose to use the assumption that the rates of information propagation along the pathway is homogeneous and, instead, letting the length of the pathway be a free parameter. We first focus on linear pathways that are sequential and have first-order kinetics, and we show how this assumption results in a three-parameter model that consistently outperforms its truncated rival and a delay differential equation alternative in recapitulating observed dynamics. We then show how the proposed assumption allows for similarly terse and effective models of non-linear pathways. Our results provide a foundation for well-informed decision making during model simplifications.

Auxin transport model for leaf venation.

The plant hormone auxin controls many aspects of the development of plants. One striking dynamical feature is the self-organization of leaf venation patterns which is driven by high levels of auxin within vein cells. The auxin transport is mediated by specialized membrane-localized proteins. Many venation models have been based on polarly localized efflux-mediator proteins of the PIN family. Here, we investigate a modelling framework for auxin transport with a positive feedback between auxin fluxes and transport capacities that are not necessarily polar, i.e. directional across a cell wall. Our approach is derived from a discrete graph-based model for biological transportation networks, where cells are represented by graph nodes and intercellular membranes by edges. The edges are not a priori oriented and the direction of auxin flow is determined by its concentration gradient along the edge. We prove global existence of solutions to the model and the validity of Murray's Law for its steady states. Moreover, we demonstrate with numerical simulations that the model is able connect an auxin source-sink pair with a mid-vein and that it can also produce branching vein patterns. A significant innovative aspect of our approach is that it allows the passage to a formal macroscopic limit which can be extended to include network growth. We perform mathematical analysis of the macroscopic formulation, showing the global existence of weak solutions for an appropriate parameter range.

Quantitative analysis of auxin sensing in leaf primordia argues against proposed role in regulating leaf dorsoventrality.

Dorsoventrality in leaves has been shown to depend on the pre-patterned expression of KANADI and HD-ZIPIII genes within the plant shoot apical meristem (SAM). However, it has also been proposed that asymmetric auxin levels within initiating leaves help establish leaf polarity, based in part on observations of the DII auxin sensor. By analyzing and quantifying the expression of the R2D2 auxin sensor, we find that there is no obvious asymmetry in auxin levels during Arabidopsis leaf development. We further show that the mDII control sensor also exhibits an asymmetry in expression in developing leaf primordia early on, while it becomes more symmetric at a later developmental stage as reported previously. Together with other recent findings, our results argue against the importance of auxin asymmetry in establishing leaf polarity.

Crack nucleation and propagation in microcrystalline-cellulose based granules subject to uniaxial and triaxial load.

Cracking patterns in four kinds of granules, based on the common pharmaceutical excipient microcrystalline cellulose (MCC) and subject to compressive load, were examined. The initial pore structure and the location of initial failure under uniaxial compression were assessed using X-ray micro-computed tomography, whereas contact force development and onset of cracking under more complex compressive load were examined using a triaxial testing apparatus. Smoothed particle hydrodynamics (SPH) simulations were employed for numerical analysis of the stress distributions prior to cracking. For granules subject to uniaxial compression, initial cracking always occurred along the meridian and the precise location of the crack depended on the pore structure. Likewise, for granules subject to triaxial compression, the fracture plane of the primary crack was generally parallel to the dominant loading direction. The occurrence of cracking was highly dependent on the triaxiality ratio, i.e. the ratio between the punch displacements in the secondary and dominant loading directions. Compressive stresses in the lateral directions, induced by triaxial compression, prevented crack opening and fragmentation of the granule, something that could be verified by simulations. These results provide corroboration as well as further insights into previously observed differences between confined and unconfined compression of granular media.

Anisotropic growth is achieved through the additive mechanical effect of material anisotropy and elastic asymmetry.

Fast directional growth is a necessity for the young seedling; after germination, it needs to quickly penetrate the soil to begin its autotrophic life. In most dicot plants, this rapid escape is due to the anisotropic elongation of the hypocotyl, the columnar organ between the root and the shoot meristems. Anisotropic growth is common in plant organs and is canonically attributed to cell wall anisotropy produced by oriented cellulose fibers. Recently, a mechanism based on asymmetric pectin-based cell wall elasticity has been proposed. Here we present a harmonizing model for anisotropic growth control in the dark-grown Arabidopsis thaliana hypocotyl: basic anisotropic information is provided by cellulose orientation) and additive anisotropic information is provided by pectin-based elastic asymmetry in the epidermis. We quantitatively show that hypocotyl elongation is anisotropic starting at germination. We present experimental evidence for pectin biochemical differences and wall mechanics providing important growth regulation in the hypocotyl. Lastly, our in silico modelling experiments indicate an additive collaboration between pectin biochemistry and cellulose orientation in promoting anisotropic growth.

The interaction of transcription factors controls the spatial layout of plant aerial stem cell niches.

The plant shoot apical meristem holds a stem cell niche from which all aerial organs originate. Using a computational approach we show that a mixture of monomers and heterodimers of the transcription factors WUSCHEL and HAIRY MERISTEM is sufficient to pattern the stem cell niche, and predict that immobile heterodimers form a regulatory "pocket" surrounding the stem cells. The model achieves to reproduce an array of perturbations, including mutants and tissue size modifications. We also show its ability to reproduce the recently observed dynamical shift of the stem cell niche during the development of an axillary meristem. The work integrates recent experimental results to answer the longstanding question of how the asymmetry of expression between the stem cell marker CLAVATA3 and its activator WUSCHEL is achieved, and recent findings of plasticity in the system.